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ABSTRACT 

Context. The composition of Galactic Cosmic Rays (GCR) presents strong similarities to the standard (cosmic) com- 
position, but also noticeable differences, the most important being the high isotopic ratio of ^^Ne/^'^Ne, which is ^^5 
times higher in GCR than in the Sun. This ratio provides key information on the GCR origin. 

Aims. We investigate the idea that GCR are accelerated by the forward shocks of supernova explosions, as they run 
through the presupernova winds of the massive stars and through the interstellar medium. 

Methods. We use detailed wind and core yields of rotating and non-rotating models of massive stars with mass loss, as 
well as simple models for the properties of the forward shock and of the circumstellar medium. 

Results. We find that the observed GCR '^■^Ne/'^^Ne ratio can be explained if GCR are accelerated only during the 
early Sedov phase, for shock velocities >1500-1900 km/s. The acceleration efficiency is found to be of the order of 
10~®-10~^, i.e. a few particles out of a million encountered by the shock escape the SN at GCR energies. We also show 
quantitatively that the widely publicized idea that GCR are accelerated in superbubbles fails to account for the high 
2^Ne/^"Ne ratio in GCR. 
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Supernova (SN) shocks are generally thought to be the main 
accelerator of the bulk of Galactic Cosmic Rays (GCR). 
Indeed, the power of GCR in the Milky Way is estimated to 
10^^ ergs/s, corresponding to 10-20 % of the kinetic power 
of Galactic supernovae (assuming canonical values of 2 SN 
per century, each one releasing 10^^ ergs of kinetic energy). 

The site of the acceleration of GCR remains debatable 
today, despite more than five decades of theoretical and 
observational studies (e.g. Strong et al. 2007 and references 
therein). Over the years, it has been suggested that GCR 
are accelerated in 1) SN remnants (either by the forward 
or the reverse shock or both), 2) the interstellar medium 
(ISM), 3) the winds of massive stars, 4) the interiors of 
superbubbles, excavated by the massive star winds and the 
subsequent SN explosions of an OB association. 

Each one of the proposed sites has its own advantages 
and shortcomings, regarding the energetics and/or the com- 
' position of accelerated matter. For instance, it has been ar- 
, gued that the hot, low density environment of a superbub- 
ble minmizes radiative losses of SN shocks and energy losses 
of accelerated particles, thus allowing the latter to reach 
substantial energies, up to the "knee" of the GCR spectrum 
(e.g. Parizot et al. 2004). On the other hand, reverse SN 
shocks running into the SN interior carry insufficient en- 
ergy to explain the bulk GCR energetics (Ramaty et al. 
1997). Moreover, they should accelerate ^^Ni, a product of 
explosive nucleosynthesis which is unstable to e~-capture 
(with a lifetime of 10"* yr) and which has not been detected 
in GCR (Wiedenbeck et al. 1999), while a rapid accelera- 
tion would render it practically stable and thus detectable. 

It was realized early on that the elemental composition 
of GCR differs significantly from the one of the ISM. Those 



diff'erences may provide valuable information on the ori- 
gin of GCR particles and, perhaps, on the site - and even 
the mechanism - of acceleration. Volatiles behave differently 
from refractories: the former display a mass-dependent en- 
richment with respect to H, which reaches a factor of 10 for 
the heaviest of them; the latter are all overabundant (w.r.t. 
H) by a factor of 20, while C and O display intermediate 
overabundances, by factors of 9 and 5, respectively (e.g. 
Wiedenbeck 2007 and references therein). 

This complex pattern is now thought to result not from 
ionization effects (as suggested in Casse and Goret 1978, 
and further developped by Meyer 1985) but rather from ef- 
fects related to elemental condensation temperature (Meyer 
et al. 1997): refractories are locked in dust grains, which are 
sputtered by repeated SN shocks and the released ions are 
easily picked- up and accelerated (Ellison et al. 1997). This, 
quite elaborate, scheme, which builds on earlier ideas by 
e.g. Cesarsky and Bibring (1981), accounts quantitatively 
for most of the observed features of GCR source compo- 
sition; still, it leaves unanswered the key issue about the 
acceleration site of GCR (and how it affects the composi- 
tion of accelerated matter) . 

The most conspicuous feature of GCR source compo- 
sition is undoubtely the high isotopic ^^Ne/^°Ne ratio. Its 
value was measured since the late 1970ies (Garcia-Munoz, 
Simpson and Wefel 1979, Wiedenbeck and Greiner 1981). 
The most accurate measurement today, obtained from anal- 
ysis of the CRIS instrument, leads to a best estimate (Binns 
et al. 2008) of 0.387 ± 0.007 (statistical) ± 0.022 (system- 
atic). This is 5.3±0.3 times the value of the (^^Ne/^^Ne)© 
ratio in the solar wind. Contrary to the case of the ele- 
mental source GCR abundances, which may be affected by 
various physico-chemical factors (first ionization potential, 
condensation temperature, etc.) isotopic ratios can only be 
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affected by nucleosynthetic processes and thus provide cru- 
cial information on the origin of cosmic ray particles. It 
should be noticed that up to now there is now clear evidence 
for any other GCR isotopic ratio to differ from solar, with 
the potential exception of ^^Fe/^^Fe, which is estimated to 
1.5±0.3 times solar (Binns et al. 2008). 

Soon after the discovery of the anomalous GCR 
22Ne/20Ne ratio, Casse and Paul (1982) suggested that 
it could be explained by a mixture of ~2% of mate- 
rial from the wind of a WC star to 98% of material 
with standard composition. In early He-burning, ^^N (pro- 
duced through the CNO cycle in the previous H-burning 
phase) is transformed almost totally in ^^Ne through 
"N(a,7)i^F(^+)i«0(a,7)22Ne. He-burning products (like 
""^^C and ^^Ne) are expelled by the stellar winds of massive 
stars during their WC phase. The observed GCR ^^Ne/^^Ne 
ratio is obtained by assuming dilution of WC material with 
matter of standard composition. Subsequent studies put 
the aforementionned idea on a quantitative basis, with the 
use of detailed models of the evolution and nucleosynthesis 
of massive, mass losing stars (Maeder 1983, Prantzos 1984, 
Meyer 1985, Prantzos et al. 1987). In those studies, the ac- 
celeration site of GCR was considered as decoupled from 
the nucleosynthesis site, and unrelated to the fraction of 
admixtured WC material. 

Higdon and Lingenfelter (2003) evaluated quantitatively 
the ^^Ne/^^Ne ratio within a superbubble, created by the 
collective action of stellar winds and SN Shockwaves. They 
adopted stellar wind yields for ^^Ne and ^^Ne from the mod- 
els of Schaller et al. (1992) and SN yields from the mod- 
els of Woosley and Weaver (1995). They found that the 
^^Ne/^^Ne ratio in the superbubble decreases with time 
(since ^^Ne from the winds dominates the evolution of 
^^Ne/^^Ne at early times) and that its time average value 
is compatible with the GCR source ^^Ne/^^Ne inferred 
from observations. In a subsequent paper, Lingenfelter and 
Higdon (2007) recognised that the Schaller et al. (1992) 
yields of ^^Ne were highly overestimatecQ and, conse- 
quently, "... new detailed calculations of the expected GCR 
isotopic ratio are called for...", but they did not attempt 
such a re-evaluation. In the meantime, Binns et al. (2005), 
using updated wind yields of massive stars with rotation 
(from the Geneva group, see Sec. 2.2), found good agree- 
ment between the observed ^^Ne/^^Ne ratio and an admix- 
ture of ~20% material from WR stars with 80% material 
of standard composition. According to Binns et al. (2008), 
since WR stars are evolutionary products of OB stars, such 
an agreement "...suggests that OB associations within su- 
perbubbles are the likely source of at least a substantial 
fraction of GCR" . 

Howewer, theoretical studies in the past 10 
years are based mostly on the paradigm of GCR 
being accelerated in SN remnants, not in su- 
perbubbles. 
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Ptuski n fc Zirakashvili (2005^ 
Berezhko et al. 2009' 



Bercz h ko fc V51k (2006)1 , 

Huskin et ai. (2010)[ [Caprioh, Amato fc Blasi (2010 



Schure et ai. (2010)[ lEllison fc Bykov (2011)| and refer 
ences therein. 'I'he kinetic energy of the bulk motion of 
the forward shock of the SN explosion is converted to 
GCR energy through diffusive shock acceleration. The 
process is highly non-linear and involves the dynamical 



^ The reason was the excessively high mass loss rates adopted 
in that work. 



reaction of both the accelerated particles and of the 
magnetic field on the system. Those studies usually take 
into account the fact that the SN explosion often occurs 
within the cavity excavated in the interstellar medium 
(ISM) by the wind of the mas sive star prior to the 
explosion (Biermann et al. (2001)1; however, the structure 
of the circumstellar environment in that case is quite 
complex and simplified models are used for its description. 
Although Caprioli et al. (2011) considered the composition 
of GCR (H,He, CNO, MgHiAl, Fe) resulting from such 
an acceleration site, none of those studies considered the 
22Ne/20Ne ratio. 



In this work we study the ^^Ne/^^Ne ratio of GCR ac- 
celerated by the forward shocks of SN explosions, as they 
run through the presupernova winds of massive stars and 
through the interstellar medium. We consider the whole 
mass spectrum of massive stars (from ^10 to 120 Mq), in- 
cluding stars with either small or large mass losses prior 
to their explosions. We consider stellar properties (masses 
of winds, ejecta, yields etc.) from recent models with mass 
loss and or without rotation (from Hirschi et al. 2005 and 
Limongi and ChiefS 2006, respectively), the former having 
larger ^^Ne enhancements in their winds. We adopt a sim- 
plified prescription (suggested in Ptuskin and Zirakashvili 
2005 and reformulated in Caprioli 2011) to describe the 
structure of the circumstellar medium at the time of the 
explosion and we consider that GCR start being acceler- 
ated in the Sedov- Taylor (ST) phase of the SN remnant 
(see e.g. Ptuskin et al. 2010). By requiring the result- 
ing IMF averaged ^^Ne/^^Ne ratio to equal the observed 
one i?of,s = (^^Ne/2ONe)GCfl/(^^Ne/2ONe)0=5.3±O.3 we are 
able to constrain the forward shock velocity to values >1900 
km/s for rotating stars (and to >2400 km/s for non rotat- 
ing ones), i.e. we find that GCR are accelerated during the 
early ST phase, lasting for a few 10^ yr. Assuming that 
10% of the SN kinetic energy is converted to GCR, we find 
that during the acceleration period a few particles out of a 
million encountered by the forward shock are accelerated. 
Finally, we reassess the superbubble paradigm for the ori- 
gin of GCR, by evaluating consistently the ^^Ne/^'^Ne ratio 
with the aforementioned stellar yields. We find that it can 
not be as high as observed, unless some extremely favor- 
able assumptions are made (only the early period of the 
superbubble lifetime considered, no gas left over from the 
formation of the OB association). We conclude that super- 
bubbles cannot be at the origin of the bulk of GCR. 



The plan of the paper is as follows. In Sec. 2 we present 
the general " set-up" of our model: the adopted stellar mod- 
els (Sec. 2.2) and their wind yields (Sec. 2.3), the descrip- 
tion of the circumstellar environment (Sec. 2.4) and the 
evolution of the forward shock in the ST phase (Sec. 2.5). 
In Sec. 3 we present our results for the (time-dependent) 
composition of the accelerated particles, the limits imposed 
on the shock velocity by the observed ^^Ne/^^Ne ratio and 
the efficiency of the particle acceleration. Finally, in Sec. 
4 we explore the ^^Ne/^^Ne ratio of GCR, assumed to be 
accelerated inside a superbubble, and we show that it can- 
not match the oberved one (unless extreme assumptions are 
made). The results are summarized in Sec. 5. 
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2. A toy model for the composition of CR 
accelerated in massive star winds 

2.1. The set-up 

The method adopted here in order to calculate the com- 
position of matter accelerated by a single SN explosion is 
schematically illustrated in Fig. [TJ At the end of its life and 
at the time of its SN explosion, a star of initial mass M* is 
left with a mass Mexp, surrounded by a circumstellar shell 
of mass Mwind=M^*-MExp, which has been lost through 
stellar wind during its prior hydrostatic evolution After the 
SN explosion, a mass of ejecta MEj=MExp-MRem (where 
Mjiem is the mass of the compact remmant, neutron star 
or black hole) expands first within the shell of mass Mwind 
and then in the ISM, with the forward shock having initial 
velocity vq = y/2Eo/MEj, where Eq is the kinetic energy 
of the SN explosion. 

In the case of stars ending their lives as WR stars, the 
wind contains both the original (nuclearly unprocessed) 
envelope of mass Meuv, and nuclearly processed layers 
of mass Mproc, enriched in products of H-burning, and 
in some cases of He-burning as well. For those stars, 
Mwind=Mproc+MEnv and Mproc is Calculated as the dif- 
ference between the mass of the nuclearly processed core 
MhcC and the mass at the explosion: Mproc=MHeC-MExp- 
For lower mass stars, exploding as red supergiants, the wind 
composition results essentially from the 1st dredge-up, i.e. 
it is a mixture of H-burning products from the stellar core 
with the original envelope composition (i.e. mass loss has 
not uncovered the He-core at the time of the explosion). 
The limit between the two classes of stars depends on their 
initial mass, mass loss rate and rotational velocity and it is 
rather poorly known at present: in general, in models with 
no rotation stars with M*>32-35 Mq become WR stars e.g. 
(Heger et al. 2002), while in models with rotation that limit 
may be as low as 22 Mq (Meynet and Maeder 2000). 

The first phase of the supernova remnant ("free expan- 
sion") takes place at shock velocity v ^ const, and ends 
when a mass Msi^Mej has been swept up in front of the 
shock wave, at which point the ST phase sets in. Following 
Ptuskin et al. (2010), we assume that efficient GCR accel- 
eration starts at this time, where the situation is energeti- 
cally most favorable. In our baseline model we shall consider 
constant acceleration efficiency ; time-dependent efficiency 
of particle acceleration is the subject of current researches 
(see Ellison and Bykov 2011, Drury 2011, and references 
therein) and will be briefly discussed in Sec. 3.3. 

The ST phase proceeds adiabatically, i.e. at ^constant 
energy and with decreasing velocity, until the temperature 
of the gas engulfed by the shock front drops to levels al- 
lowing a signiflcant fraction (about 50%) of the remaining 
energy to be radiated away. At that time, an amount of 
matter Ms2>>MEj has been swept up and the shock en- 
ters the "snow- plow" phase. At this point - and, perhaps, 
even earlier, during the ST phase - the forward shock is too 
weak to accelerate particles to GCR energies any more. 

In the aforementioned scenario, GCR are accelerated 
from a pool of particles with composition characteristic of 
the mass Mwind early on. Depending on the initial stellar 
mass, this composition may be rich in products of H- (and 
He-) burning. It is progressively diluted with ambient (first 
wind - with normal ^^Ne- and then interstellar) gas and at 
the end of the ST phase it ressembles closely the one of 




Fig. 1. Schematic representation of a supernova exploding 
in the wind of its parent star. The star, of initial mass 
M, explodes with a mass Mexp, i-e. it has lost a mass 
M.yjind=M^-MExp- The most massive stars become WR 
stars and their wind expels not only the H-envelope (of 
mass Metiv, with composition similar to the one of the 
ISM) but also nuclearly processed layers, of mass Mproc, 
i.e. Mwind^MEnv+Mproc, whcrc Mproc=MHeC-MExp and 
MhcC is the mass of the (H-exhausted) He-core. The star 
leaves a remnant (neutron star or black hole) of mass Mpem] 
the mass ejected in the SN explosion is MEj=MExp-Miiem- 
Efficient GCR acceleration presumably starts at the begin- 
ning of the ST phase, when a mass M51 ■^M^;^ is swept 
up in front of the SN shock wave (which is indicated by 
arrows) . 



the ISM. The GCR source composition observed on Earth 
should correspond to the average composition between the 
early ST phase and some later evolutionary stage of the 
remnant, and should result from the whole mass spectrum 
of exploding stars, i.e. it should be averaged over a stellar 
initial mass function (IMF). 

2.2. Properties of mass losing stars 

We adopt two sets of stellar models in this work. They are 
calculated for stars of solar metallicity, and in both cases 
the solar mixture of Anders and Grevesse (1979) is adopted. 
The corresponding metallicity is Z0=O.O19, substantially 
larger than more recent values (Lodders 2003, Asplund et 
al. 2010) and this difference results in particular from the 
reduction in the past decade of the solar abundances of C, 
N, O and Ne, which are key elements for the purpose of 
this work. For obvious consistency reasons, we keep here 
the Anders and Grevesse (1979) values, when comparing 
our results for GCR to solar ones. 

The first set of stellar models is the one of the Frascati 
group (Limongi and Chieffi 2006, hereafter LC06). It con- 
cerns 15 model stars between 11 and 120 Mq with mass 
loss but no rotation. The model includes all stages of hy- 
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Fig. 2. Input data adopted in this work, for solar metallic- 
ity massive stars with mass loss.Le/t: non rotating models 
from Limongi and ChiefR (2006); Right: Rotating models 
from Hirschi, Meynet and Maeder (2005). Top: Mass of 
the (H-exhausted) He core MhcCi of the star at explosion 
Mexp and of the compact remnant Mnem as a function of 
the initial stellar mass M*. Bottom: derived quantities are 
the mass of the wind Mwind=M^,-MExp, the mass of the 
nuclearly processed layer Mproc=MHeC-Miiem and of the 
unprocessed (^^Ne-normal) envelope MEnv=M^wind-Mproc 
and the mass of the SN ejecta MEj=MExp-Mjiem In all 
cases, curves are smoothly interpolated between model re- 
sults. 



drostatic nuclear burning and simulates the final stellar ex- 
plosion by imparting an initial velocity to a mass coordinate 
of 1 Mq (i.e. well inside the Fe core of the stars); the mass 
cut (the limit separating the ejecta of mass Mej from the 
compact remnant of mass Mnem) is chosen such as 0.1 
of ^^Ni is ejected by the explosion. 

The various masses involved in the "toy model" of 
Sec. 12.11 and Fig. [1] are provided in Table 2 of LC06 and 
are displayed in Fig. [2] {top left) of this work, whereas de- 
rived quantities are displayed in Fig. [5] ( bottom left) . It can 
be seen that stars with M*<20 have lost a negligible 
amount of mass prior to the explosion {Mwind<MEj) and 
the ST phase starts within the ambient ISM. Stars with 
M*>30 M0 have MHeC>MExp, i-e. they explode as WR 
stars, having expelled nuclearly processed layers in their 
winds. But only for the most massive stars (M*>60 Mq) 
one has MEj<Mproc, i-e. in the beginning of the ST phase 



the shock wave still expands into material with composition 
reflecting the one of the nuclearly processed core; for lower 
mass stars (3O<Af*/M0<6O) the ST phase starts when the 
shock wave encounters material with mixed composition 
from the core and the envelope, i.e. less enriched in H- and 
He- burning products. An interesting feature of those mod- 
els is that the mass of the ejecta AIej^W Mq, is similar for 
all of them, leading to similar properties(duration, swept- 
up mass) for their corresponding ST phases. 

The second set of stellar models is the one of the Geneva 
group, calculated by Hirschi, Meynet and Maeder (2005, 
hereafter HMM05) and includes both mass loss and rota- 
tion for stars with mass between 12 and 60 Mq; it has 
been complemented with results for stars of 85 and 120 
Mq, kindly provided by G. Meynet (private communica- 
tion). The initial rotational velocity is Uijot =300 km/s on 
the ZAMS, corrrsponding to an average velocity of 220 
km/s on the main sequence, i.e. close to the average ob- 
served value. Non-rotating models with the same physical 
ingredients (for convection, mass loss, etc.) have also been 
calculated, for comparison. The evolution has been calcu- 
lated to the end of Si-burning but the final SN explosion 
was not considered. Instead, an empirical prescription was 
used to evaluate the mass of the compact remnant. The 
masses of our "toy model", as given in Table 2 of HMM05, 
are displayed in Fig. [21 There is an important difference 
with respect to Table 2 of HMM05: they provide the mass 
of He-core at the time of the explosion and for the most 
massive stars this coincides with Mexp (i-e. the mass left 
to the star at explosion is smaller than the maximum ex- 
tent of the He-core); however, we are interested at the true 
value of MhcC (since this will determine how much mass of 
processed material the shock wave will encounter) and that 
value is obtained through the detailed results of HMM05 
(displaying the wind composition as a function of time - or 
of mass left). 

Comparing the results of LC06 and HMM05 one sees 
that rotation increases the mass loss {Mwind larger in 
HMM05), thus leaving the star with a smaller mass at ex- 
plosion {Mexp smaller in HMM05). Rotation also increases 
the size of nuclearly processed regions {MhcC larger in 
HMM05), since matter is rotationally mixed outwards to 
larger distances than achieved through convection. In turn, 
this leads to larger amounts of processed material Mproc 
for the HMM05 models. 

The aforementioned features of rotating vs non-rotating 
models, which are explained in detail in e.g. Maeder and 
Meynet (2000) are crucial in understanding the differences 
in the corresponding wind yields of the stars. 

2.3. The wind composition of massive stars 

LC06 provided (private communication) yields yi of all sta- 
ble nuclear species, from H to Ge, included in their models 
and ejected through the winds of the stars, up to the mo- 
ment of the explosion. HMM05 provide (Table 3 in their 
paper) the net yields yn.i of the winds of their models for 
stable species from ■^He to ^■^Na, from which the yields can 
be recovered through 

yi{M,) = ynAM*) + Mh^^^Xq,; (1) 

where the adopted solar values Xq.; are displayed in Table 
1 of HMM05. 
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Fig. 3. Total masses of selected isotopes in the winds of the 
massive stars, from the two sets of models adopted in this 
work. Solid curves are yields from HMM05 (with rotation) 
and dotted curves from LC06 (no rotation), in both cases 
interpolated between model results. The points correspond 
to results of the HMM05 models with no rotation and, in 
general, agree well with LC06 results except for the 60 M0 
star. 

The wind yields of a few selected species appear in 
Fig. [21 for the non- rotating models of LC06 and for both the 
non-rotating and the rotating models of HMM05. It can be 
seen that, in general, there is excellent agreement between 
the results for non-rotating models of LC06 and HMM05, 
for stars up to 40 Mq. Their results differ only for the 60 
M© model (and presumably for higher masses as well) and 
only for the cases of the He-burning products ^^C, ^^0 and 
^^Ne. Since both HMM05 and LC06 use the same prescrip- 
tions for mass loss, the reason of that discrepancy could be 
the use of a small amount of overshooting in the case of 
HMM05. 

Rotation has a twofold effect on stellar yields: it in- 
creases the size of the nuclearly processed layers (since it 
mixes material further than convection alone) and reduces 
the escape velocity in the stellar equator, allowing larger 
amounts of mass to be ejected in the wind. Both effects 
enhance the wind yields up to some mass limit; above it, 
the wind has removed so much mass, that less material is 
left in the star to be processed in subsequent stages of the 
evolution, thus reducing the corresponding yields. This is 
the case, for instance, with the He-burning products ^^C 
and ^^Ne, the yields of which decrease above ^60 Mq in 



the rotating HMM05 models (see Fig. [3] and HMM05 for 
details). 

In the following we assume that the wind interaction 
with the ISM has not substantially changed the wind strat- 
ification: the forward shock will first encounter the inner- 
most wind layers, containing processed material in the case 
of the most massive stars; later it will encounter the outer 
wind layers (containing mostly the initial composition) , be- 
fore running into the ISM. Fig. 4 displays the mass inte- 
grated composition of the wind (for a few key metals), as 
encountered by the forward shock, moving outwards from 
Mexp, in the case of two rotating model stars of 25 M© and 
60 Mq (from HMM05). The quantity 



m,{M) = / Xw^ndAM) dM 

J Mexp 



(2) 



is displayed as a function of mass coordinate M, y^wind,i{M) 
being the mass fraction of isotope i in the wind of 
star of mass M*. Obviously, one has mi{MExp)—^- and 
mi{M<t)=yi{M^) , i.e. at crossing the last (outermost) wind 
layer, the forward shock has encountered the totality of the 
yield y^{M^). 

An inspection of Fig. 2] shows that: 

- in the case of the 25 M© star (no He-burning products 
encountered by the shock wave), the innermost layers con- 
tain most of the ^''N produced by the CNO cycle (its mass 
increases slowly in the outermost layers), and little of ^^C 
and ^^O, which are depleted by the CNO cycle (their mass 
increases rapidly in the outermost layers); ^°Ne and ^^Ne 
are little affected by H-burning and their integrated wind 
mass increases in a way intermediate between the cases of 

and i^C. 

- in the case of the 60 M© star (He-burning, then 
H-burning products encountered by the shock wave), the 
quasi-totality of the "'^^C and ^^Ne produced by He-burning 
are encountered in the inner layers (inside ^25 M©), while 
essentially no "'^'^N is left in that region; the majority of 
^^N is found in the region 25< M/M©<50 while substan- 
tial ^®0 is found in the unprocessed envelope (beyond 50 
M©). Finally, ^°Ne is little affected everywhere (■^const. 
mass fraction) and its integrated mass rises regularly with 
the wind mass of the star. 

If detailed information (i.e. wind composition as func- 
tion of wind mass coordinate) is available, the composition 
of the material swept-up by the forward shock can be cal- 
culated in a consistent way. In case only the wind yields t/i 
and the various characteristic masses of Fig. 2 are available, 
one may adopt the following approximations: 

(i) Stars never reaching the WR stage (Af,<22 M© for 
rotating models and Af,<34 M©for non rotating ones) have 
their envelopes fully mixed after the 1st dredge-up episode 
in the red supergiant phase; the wind composition is then 
simply obtained as: 



X,, 



(3) 



(ii) Stars reaching the WR stage eject first their 
unprocessed envelope of mass Meuv and composition 
XeTii),i=X/5'jv/,i and then the innermost processed layers of 
mass Mproc and composition ^pRoc,i such that 



AM.) 



Xf 



^Mproc{M*) (4) 
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Fig. 4. Composition encountered by the forward shock 
wave as it moves (from left to right) through the winds of 25 
M© (left) and 60 Mq (right) rotating stars. The shock wave 
starts at Mexp=W MQ&nd 14.8 Mq, respectively, and the 
ST phase sets in when a mass Msi=Mexp has been swept 
up and, presumably, well mixed. The quantities displayed 
correspond to the mass integrated composition of the wind, 
starting at Mexp- 



which allows one to derive l^prod, i-e. an average abun- 
dance in the processed layers, and obtain thus an average 
wind composition profile 

XwindA^) = XprocAM*) forMExp <M < MhcC (5) 
XwindAM) = XEnv.iiM*) forMffeC < M < M* (6) 

The method outlined here turns out to provide a good 
approximation for the detailed composition profiles of the 
HMM05 models and has been used here for both sets of 
LC06 and HMM05 models, for consistency reasons. Taking 
ito account all the uncertainties of the models (prescrip- 
tions for mass loss, convective and rotational mixing etc.), 
we consider the aforementioned approximation as fairly sat- 
isfactory for the purpose of this work. 

2.4. The circumstellar environment of massive stars 

According to the ideas outlined in Sec. 2.1 particle accel- 
eration to GCR energies starts when a mass Msi ^Mej 
has been swept-up by the forward shock. The composition 
of that material, at the very beginning of particle accelera- 
tion, is the one of the material located at mass coordinate 
Ai=Mexp+ Mej and constitutes obviously an extreme for 
the abundances of the corresponding elements (upper limit 
for those produced in the stellar interior and lower limit for 
those destroyed, like H). In the following we shall assume 
that at the time of the explosion the stellar wind has fully 
kept the stratification of its various layers, as they were 
progressively leaving the stellar surface. 

The circumstellar region hit by the forward shock has 
a complex structure, depending on the properties of the 
exploded star. That structure has been explored in some 
detail with hydrodynamical models by Garcia-Segura et al. 
(1996a, b) for non-rotating stars. In the case study of a 35 



M© star, they find that the rapid wind of the O star exca- 
vates a large bubble (radius ^36 pc) of low density (10""^ 
cm""^). Inside it propagates a dense (~1 cm~'^), slow wind 
- released during the red supergiant phase - which occupies 
the innermost few pc (depending on its assumed velocity, of 
the order of a few tens of km/s). The subsequent fast (~10^ 
km/s) WR wind compresses the RSG wind, and most of the 
mass of the latter is found within a thin shell. 

The aforementioned study illustrates the complexity of 
the situation, but its results can hardly be generalized to 
the whole mass range of massive stars (for instance, lower 
mass stars will not display the fast WR wind). In fact, its 
results cannot even be safely used for the 35 M© star, since 
they depend so critically on the adopted parameters of the 
model (mass loss rates and wind velocities for the various 
stages). And they certainly fail to describe the situation 
for rotating stars, which display slow but intense (and not 
radially symmetric) mass losses on the Main sequence. 

In view of these uncertainties, we adopt here a simplified 
prescription for the structure of the circumstellar bubble, 
assuming spherical symmetry in all cases. We assume that 
the winds have excavated a bubble of mean density 7io=0.01 
cm~'^ and, consequently of radius 



Rb 



3 Af, 



w 



4tt Pa 



1/3 



(7) 



with po=no rup, rrip being the proton mass. Inside the bub- 
ble, the density profile is p{r) oc r^^, i.e. it corresponds to 
a steady stellar wind with mass loss rate Mw and velocity 
vw, which is given by 



p{r) 



M, 



w 



4:TTVwr 



(8) 



Our choice of po automatically fixes the p(r) profile and 
corresponds to a combination of Mw and vw values. 
Obviously, one has 



p(r) dr 



(9) 



Outside Rw we assume an ISM with constant density 
P/SM=1 Trip cm~'^. Our approach is similar to Caprioli 
(2011), but we do not consider here the more complicate 
case of a WR wind overtaking a RSG wind. 



2.5. Evolution during tlie Sedov-Taylor (ST) phase 

We follow the propagation of the forward shock first 
through the wind bubble and then through the ISM with 
a simple model presented in Ptuskin and Zirakashvili 
(2005) and, in a more concise form, in Caprioli (2011, his 
eqs. 3.4 to 3.9). We start the calculation from the "free- 
expansion" (ejecta dominated) phase, where the swept-up 
mass is smaller than Msj and which can be described 
by self-similar analytical solutions. In the subsequent ST 
phase (swept up mass < Mej), the model is based on the 
"thin shell" approximation (e.g. Ostriker and McKee 1988) 
which assumes that the swept up mass is concentrated in a 
thin shell behind the shock. We solve numerically the time- 
dependent equations for continuity of mass, energy and mo- 
mentum, to recover shock radius and velocity as a function 
of time. Unlike Caprioli (2011) we assume full adiabaticity 
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in the ST phase, i.e. we do not take into account the 10-20% 
energy losses of the shock through CR acceleration, which 
would reduce the shock velocity (cx Eq-^) by less than 10%. 
The mass swept up in the ST phase inside shock radius Rs 
is 



Ms{< Rs) = Mej + 47r / p{r)r^dr 



(10) 



Rsi 



We follow the evolution all the way through the ST phase, 
which ends when a significant fraction of the energy of the 
cooling remnant is radiated away (through recombination 
emission); for a solar mixture this occurs at time 



ts2 ^ 4.4 10^ yr 



/ Eo 



2/9 



yiO^^erg^ V cm^ / 



-5/9 



(11) 



In the framework of this simple model we are able to 
calculate the composition of the material encountered (and 
presumably accelerated) by the forward shock as a function 
of time, or of the swept-up mass: indeed, the integrated 
mass of each element swept up by the forward shock is 
given by an equation similar to Eq. 2 



m,{Ms) 



Msi<Rs) 



X,{Ms) dMs 



(12) 



Me. 



where Xi{Ms)=Xwind,i for Ms <Ah (Sec. 2.3) and 
Xi{M)='XiSM,i for Ms >M,, i.e. when the shock prop- 
agates in the ISM; the upper limit in the integral is given 
by eq. (10). 

Eq. 12 allows one to link the stellar model, i.e. the abun- 
dance profiles X(M), to the evolution during the ST phase 
through Eq. 10, and to the properties of the shock wave. 
Since the mass M52 swept-up in the end of the ST phase is 
much larger than the wind mass in all cases (a few 10'^ , 
compared to ~100 MQat most) for the largest part of the 
ST phase the swept up material has ISM composition. In 
order to obtain significant deviations from the solar com- 
position, such as the observed ^^Ne/^°Ne ratio, one should 
assume that significant acceleration occurs only in the early 
ST phase, when the forward shock is stronger and its ve- 
locity higher. 

3. Results 

3.1. Composition of matter in the ST phase 

Fig. [5] (top) displays the evolution of the velocity vs and 
radius Rs of the forward shock and of the mass Ms{< Rs) 
swept up by it, for the cases of a 20 and a 60 M© rotating 
star, respectively. The density of the unperturbed ISM is 
taken to be 1 cm^^ in all cases. 

The similarity of the curves for vs, Rs and Ms{< Rs) 
for the 20 and 60 Mq stars simply reflects the self-similarity 
of the ST solution. The small differences in the early ST 
phase are due to the difference of the ejected mass Mej 
in the two stars (6 vs 10 Mq, see Fig. 2), since the other 
parameters (£^0 and uism) are the same. 

The bottom panel of Fig. [5] displays the evolution of 
the i^'N/iSQ and 22Ne/2°Ne ratios (i.e. the ratios of the 
corresponding masses of Eq. 12 for each isotope), expressed 
in units of the corresponding solar values. In both cases, the 
forward shock first encounters layers with low ^^O and high 
^"^N, resulting in a high ^"^N/^^O ratio (in the 60 Mq star. 
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Fig. 5. Evolution of the SN remnants of two rotating stars 
of initial masses 20 Mq (dotted) and 60 Mq (solid), re- 
spectively. Top: velocity of the forward shock and swept 
up mass. Middle: Shock radius and density profile before 
shock arrival (see text). Bottom: Composition of the matter 
swept up up to time t by the shock wave, for the 
and ^^Ne/^°Ne ratios. In all panels, the thick portions of 
the curves indicate the period of efficient particle accelera- 
tion, i.e. from the beginning of the ST phase and as far as 
Vs > Vmin- Thc valuc of Vmin (=1900 km/s for the rotating 
star models in the figure) is chosen so that the IMF aver- 
aged theoretical ratio of ^^Ne/^°Ne matches the observed 
one in GCR (see text, Eq. (13) and Fig.©. 



^^N is depleted from He-burning in the innermost layers, 
resulting in a slightly smaller ratio than in the 20 

Mq case) . Subsequently, in the 20 Mq star the shock moves 
rapidly through the small remaining stellar envelope (3.4 
Mq) and starts propagating in the ISM, thus decreasing 
rapidly its mass-integrated ratio. In the 60 Mq 

star, the shock runs through ~20 Mq of processed material, 
with a high value of "N/^'^O, before getting to the ISM; the 
corresponding ratio decreases more slowly than in 

the 20 Mq case. 

The evolution of ^^Ne/^°Ne is quite different in the two 
models. In the 20 Mq star, no He-burning products are 
encountered by the shock wave and the ^^Ne/^°Ne ratio 
has always its initial (solar) value. A high ^^Ne/^'^Ne ratio 
is initially encountered in the processed layers of the 60 
Mq star, which is progressively diluted as the shock moves 
outwards. 
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Fig. 6. Abundance ratios of various nuclear species in GCR 
source normalized to the corresponding solar ones, as a 
function of the initial stellar mass. In all panels, solid curves 
correspond to models of HMM05 and dotted curves to mod- 
els of LC06. Upper, thin curves are for GCR accelerated at 
the beginning of the ST phase and lower, thick curves for 
the time-average at the end of GCR acceleration. An aver- 
age over a Salpeter IMF (and accounting for the swept up 
mass in each case) produces the vertical segments to the 
right, their top point corresponding to the beginning and 
the bottom one to the end of the GCR acceleration phase, 
respectively (also indicated hy filled squares). These results 
are compared to GCR source abundance ratios as derived 
by ACE data (points at the extreme right with error bars) 
in Binns et al. (2005). The most significant, unaffected by 
FIP, volatility etc., is the one of ^^Ne/^ONe. The end of the 
GCR acceleration phase is assumed to correspond to shock 
velocities Vmin such that the time and IMF averaged the- 
oretical ratio (squares) of ^^Ne/^"Ne matches the observed 
one (see text). For the set up adopted here we find we find 
u„j„=1900 km/s for HMM models and Um.m=2400 km/s 
for LC06 models. 



3.2. Composition of accelerated particles 

The composition curves of the bottom panel of Fig. [S] give 
the time (or mass) integrated composition encountered by 
the shock wave and, in consequence, the composition of 
the particles that have been accelerated up to that time. 
We assume here that particles are accelerated to GCR en- 
ergies with the same efficiency for shock velocities higher 
than some critical value Vmin, which is the same for all 



stellar masses. We determine fmm empirically by requir- 
ing that, when averaged over a stellar Initial mass function 
(IMF) $(M,), the ratio 



M22 _ /ioMe''^^*™22(M,)$(MO 



M2, 



20 dM,m2o{M,MAQ 



Robs 



22,0 



20,0 



(13) 



where _Ro6s=5.3±0.3 is the observationally determined 
source GCR ratio of ^^Ne/^^Ne in solar units and m22(M*) 
and m2o(A'/*) are calculated from Eq. 12 for stars of mass 
and for swept up masses Ms{v > fmin)- 

We adopt here a Salpeter IMF $(Af*) oc M"^ with 
X=2.35. The results of the procedure appear in Fig. [6] for 
a few selectedabundance ratios and for the models of both 
HMM05 {solid curves) and LC06 {dotted curves). In all pan- 
els, the upper {thin) curves correspond to the composition 
accelerated at the beginning of the ST phase (maximal pos- 
sible deviations from solar composition) . It can be seen that 
stars with mass <22 M© (for HMM05) and 32 M© (for 
LC06) display no He-burning products in their accelerated 
particles. N is overabundant in lower stellar masses, due to 
the 1st dredge-up. 

The lower {thick) curves in all panels of Fig. 6 cor- 
respond to composition accelerated up to the end of the 
acceleration period which is assumed to occur for a shock 
velocity Umm- The corresponding IMF-averaged quantities 
(between 10 and 120 Mq) are displayed on the right of the 
curves: their uppermost point corresponds to the beginning 
of the ST phase and the lower one (also indicated with a 
filled square) to the end of the acceleration period, i.e. to 
i^min- The value of fmm is found to be ~ 1900 km/s for 
the HMM05 models with rotation and --20% higher (2400 
km/s) for the LC06 models without rotation. The reason 
for that difference is, of course the fact that rotating models 
have larger processed layers, requiring more dilution with 
circumstellar material. 

The top and middle panels of Fig. 6 display the corre- 
sponding ratios for -'^^C/^^0 and ^^N/^^0, respectively. In 
both cases, the IMF averaged ratios are higher than the ob- 
served ones in GCR (which are, in their turn, higher than 
solar), by factors of 1.5-2. Unlike ^^Ne/^^Ne, these are ra- 
tios of different elements, having different atomic proper- 
ties. GCR source abundances are known to be affected by 
atomic effects, e.g. First Ionization Potential (FIP), or, per- 
haps more plausibly, volatility and mass/charge ratio (see 
extensive discussion in Meyer et al. 1997). The analysis of 
such effects is beyond the scope of this study. We simply 
notice here that observations indicate that refractory ele- 
ments are relatively more abundant in GCR sources than 
volatiles. Meyer et al. (1997) attribute that to the fact that 
refractories are locked up in grains, which are sputtered 
by the shock wave and the released ions are easily picked 
up and accelerated. ^^O being more refractory than both 
^"^N and ^^C, it is expected that its abundance in GCR 
source will be enhanced by that effect, and the correspond- 
ing 12C/160 and "N/I60 ratios in GCR source will be 
lower than predicted from stellar nucleosynthesis alone. The 
fact that in Fig. 6 we obtain higher than observed ratios for 
12(-<^16q i^N/^^O is encouraging in that respect, since 
the aforementioned atomic effects would lower those ratios, 
hopefully to their observed valuefl 



^ One might think that a combined analysis of the ^^C/^^O 
and ^^N/^^O ratios of Fig. 6 could help to constrain the atomic 
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Fig. 7. Graphical presentation of the results discussed in 
Sec. 3.3. Particle acceleration starts at the beginning of the 
ST phase, located at mass coordinate Ai=M Exp+ Msj, i.e. 
when the forward shock (FS, arrows), launched at Mexp, 
has swept up a mass Msi—Mej- Acceleration stops at mass 
coordinate A2 , corresponding in the case dicussed here to a 
shock velocity of 1900 km/s. The mass sampled by the FS 
between those two regions is Macc=A.2-Ai. For rotating 
stars with mass M >30 Mq, an increasing part of Macc 
includes nuclearly processed material {shaded aerea), while 
for rotating stars with M <18 M0, Maccc contains only 
material of ISM (=solar) somposition. 



3.3. Efficiency of particle acceleration 

The material of Sec. 3.1 and 3.2 is summarized in Fig. [7] 
for the rotating models of HMM05. The forward shock, 
launched at mass coordinate Mexp, sweeps up a mass 
AIsi ^Mexp and then starts accelerating particles, at mass 
coordinate Ai=Mexp+Msi, up to point A2 (where its ve- 
locity becomes Vmin)- For rotating stars with M <15 Mq, 
Ai lies beyond the stellar surface and only ISM is accel- 
erated. In stars with 15<M*/Mq<25, Ai lies beyond the 
processed/mixed interior MhcC and the forward shock ac- 
celerates first envelope and then ISM material. For stars 
above 35 Mq, the shock first accelerates processed mate- 
rial (hatched aerea), then the - ^^Ne normal - envelope and 
then ISM. Finally, for stars with M>70 Mq, particle ac- 
celeration ends when the shock is still within the massive 
stellar envelope. 

The mass of circumstellar material from which particles 
are accelerated is Macc = A^— A\ and lies in the range 



processes shaping the GCR source abundances. However, the 
abundances of ^^C and ^^O are affected by the still uncertain 
value of the ^^C(q,7^^0 reaction rate and are unsuitable for 
such a study. 



30-40 M0 . The mass of particles that have been accelerated 
is given by 



rriAp = ^ NiAivrip 



(14) 



where Ni is the total number of nuclei of species i, Ai the 
corresponding mass number and mp the proton mass. The 
number Ni, or rather the product NiAi can be determined 
by noticing that the total energy carried by those acceler- 
ated particles is a fraction / of the kinetic energy Eq of the 
supernova 



N,A, / E Q{E) dE = f Eo 



(15) 



where Q{E) is the spectrum of accelerated particles {Ai 
appears on the left side of Eq. 15 because energies are ex- 
pressed in energy units per nucleon). The efficiency / of 
conversion of Eq to accelerated particles escaping the su- 
pernova is estimated to / ~0.1, while another 0.1-0.2 goes 
to the acceleration of particles which are finally trapped in 
the SN (Ellison and Bykov 2011). The particle spectrum is 
often described by a power law in momentum 



Q{E) 



P_ 

/3 



-exp(-^/^o) 



(16) 



where (3 — v/c is the velocity expressed as a fraction of the 
light velocity, p is the particle momentum per nucleon, the 
factor s is usually 2< s <3 (in the case of strong shocks) 
and Eq is a cut-off energy, here taken to be 3 TeV (the 
results are insensitive to much higher values). 

Assuming that the spectrum of accelerated particles es- 
caping the SN is given by Eq. 16, one finds the efficiency 
with which particles are accelerated from the shocked cir- 
cumstellar medium, through Eq. 14 to 16: 



W 



rriAP 
Macc 



(17) 



and for rotating stars it is found to be in the range 3-6 
10~^ i.e. a few particles out of a million encountered ones 
are accelerated by the forward shock to GCR energies. In 
the case of non-rotating stars, the energetics is the same, 
but the swept-up mass is smaller (by a factor of two, on 
average, in order to get the observed ^^Ne/^°Ne ratio) and 
the corresponding efficiency is W '^10^^. These estimates 
constitute only a gross average, since the efficiency of parti- 
cle acceleration may depend on several factors, not consid- 
ered here, like e.g. the density of the circumstellar medium 
or the shock velocity - through a smoothly varying func- 
tion f{vs) instead of the Heavyside function /(> v„mi)—l 
and /(< Vmin)=0 considered here - or the shock radius at 
the time of acceleration, since particles may subsequently 
suffer adiabatic cooling before escaping (e.g. Drury 2011, 
Bykov and Ellison 2011). Notice that some of those effects 
may have opposite time dependencies. For instance, parti- 
cles accelerated earlier on (at higher shock velocities and 
presumably with higher efficiencies) are expected to suffer 
more from adiabatic cooling (because they are produced 
at smaller radii) . These effects require a much more thor- 
ough investigation. As a first step towards that direction, 
we also tested the case where the efficiency of shock accel- 
eration varies with shock velocity as / oc (Drury 2011). 
In that case, material with high ^^Ne/^°Ne is efficiently ac- 
celerated in the inner layers of the most massive stars. In 
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order to obtain again the observed GCR source ^^Ne/^°Ne 
ratio, one has to dilute the mixture by allowing acceler- 
ation for shock velocities lower then the reference values 
found above: we thus find values of f mm =1600 km/s for the 
HMM05 yields and 2150 km/s for the CL06 yields; the cor- 
responding overall acceleration efficiencies increase then by 
~40% from the reference values given above, ranging from 
4-7 10-6 for HMM05 yields to 1-1.4 10"^ for CL06 yields. 
The reason for obtaining such a small difference (only a 
few hundreds of km/s) between the non- realistic reference 
case and the - perhaps more realistic - case of velocity de- 
pendent efficiency, is the adopted unified treatment, which 
considers acceleration in both low mass and high mass su- 
pernova: the former accelerate almost pure ISM with low 
^^Ne/^°Ne and the latter almost pure wind material with 
high ^^Ne/^'^Ne. Allowing for lower values of Vmin involves 
a considerably larger amount of ISM processed by low mass 
supernovae, since the lower acceleration efficiency (/ oc u^) 
is more than compensated by the n/5M47rr^c?r factor. 

As approximate as they may be, the results obtained 
here through the constrain of the observed GCR source 
^^Ne/^°Ne ratio, indicate clearly that acceleration has to 
occur for shock velocities larger than ^^1500 km/s, and they 
may help to improve our understanding of particle acceler- 
ation in SN remnants. 



4. GCR cannot be accelerated (mainly) in 
superbubbles 

The idea that CGR are accelerated mainly in superbubbles 
has been suggested by Kafatos et al. (1981) and reassessed 
by Higdon et al. (1998). Massive stars are mainly formed in 
OB associations and the winds of the most massive of them 
initiate the formation of superbubbles, while the subsequent 
SN explosions power the expansion of those superbubbles 
for a few 10'' years. Higdon et al. (1998) argued that the 
environment of such superbubbles, enriched with the ejecta 
of stellar winds and core collapse SN explosions, provides a 
composition that compares favourably to the inferred GCR 
source composition. 

One potential problem with that idea is that ^^Ni, a 
well known product of SN nucleosynthesis, is absent from 
GCR arriving on Earth (Wiedenbeck et al. 1999). Since 
^^Ni is unstable to electron capure, with a lifetime of ^10^ 
y, Higdon et al. (1998) argued that SN explosions occur 
within the superbubble with a sufficiently low frequency 
(less than one SN every 3 10^ y), as to allow for the decay 
of 59]\[i between two SN explosions. Prantzos (2005) pointed 
out that energetic stellar winds can also accelerate particles 
and they are not intermittent (like SN) but occur contin- 
uously in superbubbles; in that case ^^Ni is continuously 
accelerated to GCR energies and, being unable to capture 
an electron, it becomes effectively stable and it should be 
detectable in GCR. Binns et al. (2008) counter-argued that 
the period of energetic WR winds represents a small (al- 
beit not negligible) fraction in the early lifetime of an OB 
association and therefore only little (and presumably unde- 
tectable) amounts of ^^Ni would be accelerated, thus saving 
the "superbubble paradigm". 

In the meantime, Higdon and Lingenfelter (2003) eval- 
uated the isotopic composition of ^■^Ne/^°Ne - the critical 
abundance anomaly in GCR - expected in a superbubble, 
based on (a very heterogeneous set of) then available yields 



of WR stars and SN. They found that the GCR ^ZNe/^^Ne 
ratio "...can be easily understood as the result of GCR ac- 
celerated primarily out of superbubbles with a mean metal- 
licity ZsB=2.7±0.4 Zq..." and that this result "... pro- 
vides strong, additional evidence for a superbubble origin 
of GCR". In a subsequent paper, Linfenfelter and Higdon 
(2007) explored quantitatively that senario for other GCR 
abundances, finding again that results compare favourably 
to observations. 

Despite their apparent sophistication, the aforemen- 
tioned arguments for a superbubble origin of GCR miss a 
simple point: massive stars are the principal source of both 
^^Ne and ^^Ne in the Universe. As a result, the ^^Ne/^°Ne 
ratio of a generation of massive stars (integrated over an 
IMF and including WR winds and all kinds of SN ejecta) 
should be solaiQ. Since ^^Ne is overabundant with respect 
to ^°Ne only in the massive star winds (and not in the SN 
ejecta or the ISM), a popular exercice in the past 30 years 
or so consisted in evaluating the mixing ratio of WR ejecta 
with SN ejecta (or with average ISM), as to obtain the 
observed ^^Ne/^°Ne ratio in GCR. (see Introduction and 
references therein). 

However, it is expected that a properly weighted mix- 
ture of the ejecta of massive stars in a superbubble (i.e. in- 
cluding stellar winds and SN ejecta and folded with a stellar 
IMF) would produce a solar ^^Ne/^°Ne ratio. In this sec- 
tion, the same exercise is repeated with the yields of CL05 
and HMM05 and by taking into account the correspond- 
ing stellar lifetimes. In Fig. [5] we present again the adopted 
20Ne and ^^Ne yields of LC06 and HMM05 for the stellar 
winds (solid and dotted curves, i.e. the same as in Fig. 3) 
and also the total yields (winds plus SN ejecta, points). 
The latter are from LC06 for stars with mass loss (but no 
rotation) and from Woosley and Weaver (1995, WW95) for 
stars without mass loss; they are in good overall agreement, 
except for ^'^Ne in stars around 20 Mq. For the SN ejecta, 
we consider only stars with M<40 M0, i.e. we assume that 
more massive stars eject their ^"Ne and ^^Ne through their 
winds and then form black holes; in that case we maximize 
the ^^Ne/^'^Ne ratio expected from massive stars. 

An inspection of Fig. [5] shows that the total (SN) ^^Ne 
yields of M<40 Mq stars are comparable to the wind ^^Ne 
yields of the most massive (~100 Mq) stars; however, the 
corresponding total ^"Ne yields of M<40 Mq stars are at 
least 10 times larger than the wind ^'^Ne yields of the most 
massive (~100 Mq) stars. Thus, while the production of 
^^Ne receives a sizeable contribution from the most massive 
stellar winds (at least for the rotating stars), the production 
of 2°Ne is totally dominated by the SN ejecta of M<40 Mq 
stars. In a coeval stellar population, like the one expected 
in an OB association or a superbubble, the ^^Ne/^'^Ne ratio 
will evolve then from higher than solar to ^solar values, as 
the M<40 Mq stars eject their core products a few Myr 
after their more massive counterparts. 

We quantitative illustrate this effect in Fig. ^ where we 
present the evolution of a stellar population born at time 
t=0 with a Salpeter IMF and total mass of 1 Mq (results 
can be directly scaled to masses of OB associations, while 



^ Strictly speaking, this concerns massive stars of roughly so- 
lar metallicity, such as those that contributed to the composition 
of the solar system; however, metallicity has evolved very little 
in the past several billion years in the local volume of ~ 1-2 kpc 
radius, where most GCR originate. 
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Fig. 8. Yields of ^'^Ne (top) and ^^Ne bottom) as a function 
of the initial stellar mass. Wind yields are from HMM05 
{solid curves) and LC06 {dotted curves), respectively (same 
as in Fig. 3). Total yields are from LC06 {filled squares) and 
WW95 {asterisks). 



abundances and abundance ratios remain the same). The 
top right panel displays the rate of SN explosions; for a 
total stellar mass of 10'' M0, resulting in a ^100 OB stars, 
there are about 3 SN/Myr, i.e. -1 SN every 300 000 yr, 
as evaluated in Higdon et al. (1998) in order to avoid the 
problem of ^^Ni acceleration. The top left panel displays 
the ejection rates of ^"^Ne and ^^Ne (in Mo/Myr per Mq 
of stars formed) : ^^Ne from stellar winds dominates only 
for a couple of Myr but ~5 Myr after the stellar formation, 
^°Ne from M<40 MQexplosions dominates. Notice that, in 
order to maximise the ^^Ne/^°Ne ratio, we adopted only the 
yields of HMM05 for these calculations (i.e. corresponding 
to the solid curves in Fig. |8] with only wind yields above 
M=40 Mq). 

The instantaneous (^^Ne/^°Ne)sB ratio in the super- 
bubble at time t is the ratio of the amounts of ^°Ne and 
^^Ne cumulated up to that time, which are provided in the 
middle panels of Fig. [HI In order to evaluate those quanti- 
ties, one needs to make an assumption about the amount of 
pre-existing gas still left inside the superbubble. In the left 
middle panel of Fig. [3] it is assumed that no gas is left after 
star formation: gas is exclusively supplied by the wind and 
explosive stellar ejecta. In that case, ^^Ne dominates early 
on, but ^°Ne takes over after a couple of Myr. In the right 
middle panel it is assumed that an amount of gas with solar 
composition and equal to 20% of the mass of formed stars 
remained in the superbubble; the stellar ejecta are diluted 
into it and the superbubble composition is now dominated 
always by ^°Ne. 

The bottom panels of Fig. ^ display the corresponding 
(^^Ne/^'^Ne)sB ratios in the superbubble. In the case of 
zero initial gas (left), the instantaneous (^^Ne/^°Ne)sB ra- 
tio remains high - from 18 to 5 times solar - for a couple of 
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Fig. 9. Evolution of composition in a superbubble. Top left: 
Ejection rate of '^^Ne{dotted) and ^^Ne {solid) in Mo/Myr 
per 1 M0 of stars formed (assuming a Salpeter IMF be- 
tween 0.1 and 120 M©). Top right: Supernova rate per Myr. 
Middle: Total mass of ^°Ne {dotted) and ^^Ne {solid) accu- 
mulated in the superbubble, in case there was no initial gas 
there (Mo=0, left) and in the presence of an initial mass of 
gas equal to 1/5 of the formed stars {Mo=0. 20, right); in the 
former case ^^Ne dominates by mass early on, whereas in 
the latter the mass of ^°Ne dominates from the very begin- 
ning. Bottom: Corresponding evolution of the O/H {dashed) 
and of the {'^'^'Ne/^°'Ne)sB ratio {dotted) in the SB gas; the 
{solid curve represents the (^^Ne/^°Ne)c:i? ratio of particles 
accelerated up to that time and has to be compared to the 
observed GCR source ratio {dashed horizontal line). 

Myr, but soon after stars of M<40 Mq start dying its value 
goes rapidly to solar (as expected from the nucleosynthesis 
argument presented in the beginning of this section). The 
(^^Ne/^'^Ne)cR (< t) ratio of particles accelerated to CR 
energies up to time t is the time integral of the instanta- 
neous (^^Ne/^°Ne)sB ratio in the superbubble 

where it is assumed that the efficiency of particle acceler- 
ation does not vary with time. (^^Ne/^''Ne)cR (< t) rep- 
resents the mean CR composition accelerated by a super- 
bubble which "operates" up to time t. 

In the bottom left panel it is seen that the 
(^^Ne/^''Ne)cR (< t) ratio (in solar units) remains above 
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the observed i?obs=5.3±0.3 for about 7 Myr, while it de- 
clines steadily after that, tending to an asymptotic value of 
~2, considerably smaller than observed. The corresponding 
metallicity in the superbubble, expressed by the O /H ratio 
(dashed curve), is quite high - about 10 times solar - since 
it represents pure core collapse SN ejecta; such high metal- 
licity values have never been reported for any astrophysical 
environment (except SN remnants). 

The bottom right panel of Fig. |9] displays the corre- 
sponding quantities in the more realistic case of a superbub- 
ble endowed with some gas left over from the star forma- 
tion. In that case, although metallicity is still high (about 
3 times solar after 7 Myr), the ^^Ne/^°Ne ratios never get 
above 3 times solar, either in the superbubble or in the ac- 
celerated particles; in fact they are close to solar for the 
largest part of the superbubble lifetime. 

The results displayed in Fig. [HI suggest that the observed 
-Robs =5. 3±0. 3 ratio in GCR cannot be obtained from ma- 
terial accelerated in a superbubble: even under the most 
favorable possible conditions, such as those adopted here 
(yields from rotating massive stars, only wind yields consid- 
ered above M=40 Mq, no gas left over from star formation), 
the resulting ^^Ne/^°Ne ratio remains at high values only in 
the early evolution of the superbubble, during a small frac- 
tion of its lifetime. It is interesting to notice in that respect 
that Binns et al. (2008) suggested that no substantial par- 
ticle acceleration must occur during that early period, in or- 
der to avoid the problem of ^^Ni acceleration from the WR 
winds; thus, saving the "superbubble paradigm" for the ori- 
gin of cosmic rays, requires that acceleration occurs only in 
the late superbubble evolution, when the WR winds have 
essentially stopped and only SN inject intermittently their 
kinetic energy in the superbubble. However, we have shown 
here that during this period, the (^^Ne/^°Ne)ci? ratio is 
considerably lower than observed. In fact, under realistic 
conditions (some initial gas left, explosive yields also above 
40 Mq considered, at least for some stars), the ^^Ne/^°Ne 
ratio in a superbubble would never get to values as high as 
observed, contrary to claims made in the literature in the 
past decade. It must be stressed that this conclusion does 
not depend on detailed adopted yields, but on a simple ar- 
gument, namely that the IMF averaged ^^Ne/^°Ne ratio in 
a superbubble has to be close to solar during the largest 
period of the superbubble lifetime. 



5. Summary 

In this work we explore some implications of the idea that 
cosmic rays are accelerated by the forward shock in super- 
nova remnants, during their ST phase. We focus on the 
chemical composition resulting from such an acceleration 
and, in particular, on the ^^Ne/^^Ne ratio, which is the most 
characteristic feature of the observed GCR source composi- 
tion and is unaffected by atomic effects. For that purpose, 
we adopt recent models of the nucleosynthesis and evolu- 
tion of massive stars with mass loss: those of LC06 with no 
rotation and those of IIMM05 with rotation. 

In Sec. 2 we present a detailed summary of the proper- 
ties of those models and, in particular, of the chemical com- 
position of their winds, insisting on the fact that rotating 
models release more ■^'^Ne in their winds than non-rotating 
ones. We also present the adopted model for the evolution 
of a SN remnant within a stellar wind, based on the ideas 



of Ptuskin and Zirakashvili (2005) and the equations sum- 
marized in Caprioli (2011). 

In the framework of the simple model adopted here, 
we follow the time-dependent composition of GCR, accel- 
erated by the forward shock as it runs through either the 
ISM (in the case of stars with M<25-35 Mq, depending on 
rotation) or through the stellar wind (in the case of more 
massive stars). In fact, during the largest part of the ST 
phase, the shock runs through ISM and encounters a ~solar 
composition. In order to reproduce the observed high value 
of (22Ne/2°Ne)c_R (i?o6s=5.3±0.3 in solar units) after ac- 
counting for the stellar IMF, we have then to assume that 
acceleration is efficient only during a short early period in 
the ST phase. We chose to use the shock velocity as a crite- 
rion for efficient acceleration and, based on the aforemen- 
tioned CR composition argument, we find that shock ve- 
locities larger than ^^1900 km/s (for the rotating stellar 
models) or 2400 km/s ( for the non-rotating ones) are re- 
quired. This result is obtained by assuming a step function 
for the efficiency / of particle acceleration (/=0 before the 
ST phase and after Vmim and /=1 between the two). For 
the - perhaps, more realistic - assumption of a velocity- 
dependent efficiency f oc v'^, we find slightly lower values 
for v,nin (1600 km/s for HMM05 yields and 2150 km/s for 
CLOG yields, respectively). In the framework of the adopted 
models, this corresponds to a circumstellar mass of several 
tens of Mq encountered by the forward shock. Assuming, 
furthermore, that 10% of the SN kinetic energy is used in 
acceleration of escaping cosmic rays with standard energy 
spectra, allows us to evaluate the efficiency of that acceler- 
ation: we find that a few particles out of a million encoun- 
tered by the forward shock are accelerated to CR energies. 
We also notice that this scheme of GCR acceleration does 
not suffer from problems related to the absence of unstable 
^^Ni in observed GCR composition: this heavy nucleus is 
well inside the SN ejecta and is not reached by the forward 
shock which accelerates only wind material and ISM. 

The aforementioned scenario assumes that even the 
most massive stars, up to 120 Mq, develop strong for- 
ward shocks and accelerate the particles of their WR winds. 
For non-rotating stars, this is a rather extreme assump- 
tion, since it has been argued (Heger et al. 2003) that non- 
rotating masssive stars of about solar metallicity collapse 
into black holes, if their mass is in the 30-60 M© range (see 
their Fig. 1). Notice that stars in the 60-120 Mq range (the 
most important ^^Ne producers) end as black holes in their 
scheme. However, for slightly higher metallicities - such as 
those encountered in the inner Galactic disk - they find that 
only neutron stars are formed, because higher stellar mass 
losses result in a less massive star at explosion. It should be 
noticed that the details of massive star explosions remain 
poorly understood at present (see e.g. Hanke et al. 2011) 
and so is the fate of a massive star above 30 Mq (see Fryer 
et al. 2011 for a recent - but certainly not definitive - assess- 
ment). The fate of the rotating mass losing stars considered 
here is even less well known. 

In view of the afornientioned uncertainties, we feel that 
the scenario proposed here can be considered as valid at 
present, although future refinements in our understand- 
ing of massive star explosions may change it quantitatively 
(and even qualitatively, if it turns out that most masive 
stars above, say, 50 Mq, end up as black holes). 

Finally, we explore the idea that CR are accelerated in 
superbubbles, in which case their composition results from 



12 



N.Prantzos: On the origin and composition of Galactic cosmic rays 



the ejecta of both stellar winds and SN explosions. We first 
notice that simple nucleosynthesis arguments suggest that 
the resulting composition, averaged over the stellar IMF, 
should be very close to solar. We demonstrate this quanti- 
tatively, with a simple model for the evolving composition 
of a superbubble, enriched first by the (^^Ne rich) winds of 
the most massive stars, then by the (^°Ne rich) SN ejecta 
of less massive stars. We find that, after a few Myr the 
superbubble (^^Ne/^°Ne)5_B ratio tends to solar, and so 
does the average (^^Ne/^°Ne)cfl ratio in accelerated par- 
ticles. We conclude that superbubbles cannot provide the 
observed high Robs value of CR sources and, therefore, are 
not the main site of CR acceleration. On the contrary, SN 
remnants - including those expanding in the pre-explosion 
environment of a stellar wind - appear as suitable sites of 
GCR acceleration. 
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